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Abstract: We investigate a lattice version of the Yang-Lee model which is character- 
ized by a non-Hermitian quantum spin chain Hamiltonian. We propose a new way to 
implement T-'T-symmetry on the lattice, which serves to guarantee the reality of the spec- 
trum in certain regions of values of the coupling constants. In that region of unbroken 
"PT-symmetry we construct a Dyson map, a metric operator and find the Hermitian 
counterpart of the Hamiltonian for small values of the number of sites, both exactly and 
perturbatively. Besides the standard perturbation theory about the Hermitian part of the 



(N 
> 



Hamiltonian, we also carry out an expansion in the second coupling constant of the model. 



o 



Our constructions turns out to be unique with the sole assumption that the Dyson map is 
Hermitian. Finally we compute the magnetization of the chain in the z and x direction. 

O 

1. Introduction 

" r- i i 

^ ■ It is known for about thirty years that ordinary second order phase transitions can be 

described by the Yang-Lee model ]l], ||, [|. This model admits a quantum field theoret- 
ical description in form of a Landau-Ginzburg Hamiltonian for a scalar field 4> with an 
additional (^-interaction and a term linear in the scalar field with an imaginary coupling 
constant. The model has been identified |] as a perturbation of the A^s^-model in the 
7W Pi9 -series of minimal conformal field theories f5|. It is the simplest non-unitary model 
in this infinite class of models, which are all characterized by the condition p — q > 1 and 
whose corresponding Hamiltonians are all expected to be non-Hermitian. 

Here we shall investigate a discretised lattice version of the Yang-Lee model considered 
by von Gehlen §, 0, which is an Ising quantum spin chain in the presence of a magnetic 
field in the z-direction as well as a longitudinal imaginary field in the x-direction. The 
corresponding Hamiltonian for a chain of length N is given by 

1 - 

H(\,k) = --^(a z j + \aj(jj +1 +iKO-j), A,k£R. (1.1) 
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It acts on a Hilbert space of the form (C 2 )® N where we employed the standard notation for 
the 2 N x 2 Ar -matrices erf ,y,z = I® I ® . . . ®er x ' y,z <8> • • • <8>I <8> I with Pauli matrices describing 
spin 1/2 particles 






as ith factor acting on the site % of the chain. Their commutation relations are direct sums 
of su(2) algebras 

[a], a y k ] = 2ia z 5 jk , [a z p a%\ = 2ia y 5 jk , [a), a%] = 2icr*5 jk , with j, k = 1, . . . , N (1.3) 

A further real parameter (3 may be introduced into the model by allowing different types of 
boundary conditions <r^* = fio x ^ y ' z , albeit here we will only consider the case of periodic 
boundary conditions and take (3 = 1. 

Since all Pauli matrices are Hermitian it is obvious that H(X, k) is non-Hermitian 

H^(X,k) = H(X,-k) ^ H(X,k). (1.4) 

This poses immediately two questions: First of all, is the spectrum still real, despite the 
fact that the vital property of Hermiticity which guarantees this is given up and second 
is it still possible to formulate a meaningful quantum mechanical description associated to 
this type of Hamiltonians? These issues have attracted a considerable amount of attention 
in the last ten years, since the seminal paper by Bender and Boettcher || and meanwhile 



many satisfying answers have been found to most of them; for recent reviews sec || [l0|, 11]. 

Our manuscript is organised as follows: In section 2 we present various alternatives 
about how 'PT-symmetry can be implemented for quantum spin chains. In section 3 we 
establish our notation and recall some of the well known facts concerning a consistent 
quantum mechanical framework for 'PT-symmetric systems. We analyze the model ( |1.1| ) 
in section 4 and section 5, where the former is devoted to non-perturbative and the latter 
to perturbative results. In section 6 we compute the magnetization for the model (1.1) and 
we state our conclusions in section 7. 



2. PT-symmetry for spin chains 

Preceding the above mentioned recent activities von Gehlen found numerically || 0] that 
for certain values of the dimensional parameters A and k the eigenvalues for H (A, k) are all 
real, whereas for the remaining values they occur in complex conjugate pairs. He provided 
an easy explanation for this feature: Acting adjointly on the Hamiltonian with a spin 
rotation operator 

n = eW = U^=(I + i°*h with Sf = f>f, ^ 

has the effect of rotating the spins at each site clockwise by tt/2 in the xy-plane, such that 
the corresponding map acts as 1Z : (erf, a y , a z ) — > (—erf, erf, erf). The resulting Hamiltonian 
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is a 2^ x 2^ non-symmetric matrix with real entries given by 

2 



1 N 

h(\, K ) = hh(\, k)k~ 1 = J>f + A(T M+i - ( 2 - 2 ) 



1=1 



Its eigenvalues and those of H(\, k) are therefore either all real or occur in complex conju- 
gate pairs. This is precisely the well known behaviour one finds when H(\, k) is symmetric 
with respect to an anti-linear operator 14, 15, 16 1, which as mentioned above has 



recently attracted a lot of attention. In quantum mechanical or field theoretical models the 
anti-linear operator is commonly taken to be the PT-operator, which carries out a simul- 
taneous parity transformation P : x — ► —x and time reversal T : t — > —t. When acting on 
complex valued functions the anti-linear operator T is understood to act as complex con- 
jugation. Real eigenvalues are then found for unbroken PT-symmetry, meaning that both 
the Hamiltonian and the eigenfunctions remain invariant under PT-symmetry, whereas 
broken PT-symmetry leads to complex conjugate pairs of eigenvalues. 

We will now argue that PT-symmetry on the lattice can be interpreted in various 
ways. One may for instance reflect the chain across its midpoint via the map P' : a^' v ' z — > 
CTV+i-j as suggested by Korff and Weston |L7]] and used thereafter in [18]. It is obvious 



that the Hamiltonian (1.1) is invariant with regard to this symmetry. However, when 



keeping the interpretation of T as a complex conjugation, and thus ensuring that the P'T- 
operator is anti-linear, one easily observes that this type of transformation does not leave 



the Hamiltonian (LI) invariant, i.e. we have [P'T,H] ^ 0. 



Therefore we need to implement PT-symmetry in a different way for H(X,k) to be 



able to analyze its properties along the lines proposed in [|L2|, |13|, |b|, |15|, |lq| . We propose 
here that one carries out a parity transformation at each individual site and reflect every 
spin for instance in the xy-plane on y = —x. This is obviously achieved by 1Z 2 . As 
1Z 4 = Y\iLi(—^)i = (—1)^1®^ and not the desired identity operator, we take here 

N 

r = -in 2 = e^ sz -^=Ha!, with V 2 = I® N , (2.3) 

i=l 

as our parity operator. Consequently this transformation acts as 

V:(af,alaf)^(-al-ala!). (2.4) 
Thus with T being the usual complex conjugation, which acts on the Pauli matrices as 

T:(af,alaf)^(af,-a?,af), (2.5) 
we have identified an anti-linear operator constituting a symmetry of the Hamiltonian (|L 

[VT,H]=0. (2.6) 



This operator provides more information than the transformation ( |2.2[) , because we have 
now in addition a concrete criterium, which distinguishes the regimes of real and complex 
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eigenvalues. We can precisely separate the two domains U-p-r and Ubvi in the parameter 
space of A and k defined by the action on the eigenstates <1?(A, k) of H(X, k) 

/ $(A, K) for (A, k) G U bVT . 

According to the general reasoning provided in [l4|, |l^, 16 1, simultaneous eigenfunc- 

tions of VT and H (A, k), that is for (A, k) 6 f77?Tj are then associated with real eigenvalues 
whereas in the regime of broken "PT-symmetry, that is (A, k) £ Ubj>T, the eigenvalues 
emerge in complex conjugate pairs. 

From the above it is clear that we may define equally well different types of PT- 
operators closely related to the one introduced in ([2.3]). For instance we can define 



N N 



V x :=l[af and 7>„ : = EK' ^ 

i=i i=i 

which obviously act as 

V x :(af,ay,al)^(al-al-af) and V y : (of, of, of) - (-of , of, -of). (2.9) 



Clearly these parity operators can not be used in the same way as P in (2^3) to introduce 
a PT-symmetry for H (A, n) when keeping T unchanged. However, they serve to treat 
non-Hermitian Hamiltonians of a different kind, such as obvious modifications of H(X, k) 
and also to allow for alternative treatments of non-Hermitian spin chains, such as the 



XXZ-spin-chain in a magnetic field [19] 
N-l 



Hxxz = lY, + "Mf-i + A +(^Ui ~ !)] (2-10) 

i=l 

with A± = (q ± q~ 1 )/2 previously studied in [17, 18|. Obviously when q ^ R this Hamil- 
tonian is non-Hermitian, but we observe that it is PT-symmetric when using any of the 
parity operators defined in (|2.8|) and keeping T to be the usual complex conjugation 

[V x T,H xxz ]=0 and [V y T, H X xz] = 0. (2.11) 

Thus besides reflecting the chain across its midpoint in form of a "macro-reflections", as 
suggested in [17], we may also carry out the parity transformations on each individual side. 



It appears that these "micro-reflections" (|2.3j), (|2.8D allow for a wider range of possibilities, 



such as for instance Hamiltonians of the type H(X, k) in (1.1), which could not be tackled 
with P' : a^' y ' z — ► o"^^_^. The different possibilities are simply manifestations of the 
well known ambiguities non-Hermitian Hamiltonians possess with regard to their operator 
content This also means that the symmetries ( |2.11 ) will lead to a different kind of 



physical systems than those identified in [17] 



It is well known that Hxxz can be expressed in terms of generators of a Temperley- 
Lieb algebra Ei, i.e. simply by writing the Hamiltonian alternatively as Hxxz = X^T 1 ^i- 
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It is then trivial to see that the algebra remains invariant under a 'PT-transformation when 
realized as Q2.8Q : T : Ei — > E* , V XjV : E{ — ► E*, such that V XjV T : E{ —> Ei. On the other 
hand when implementing the "macro-reflection" on the entire chain, the 'P'T-symmetry 
on the generators is broken, i.e. V"T : E{ — > Ejf+i—i, as was found in [17]. 



A further interesting non-Hermitian quantum spin chain has recently been investigated 
by Deguchi and Ghosh [ 20 1 



N 

H DG = K zz afaf +1 + K x af + K y a\, (2.12) 

i=l 

with k zz G R and K x , n y G C. Clearly when k x or n y ^ R the Hamiltonian Hjjq is 
not Hermitian, which is the case we will consider. As the previous model also the quasi- 
Hermitian transverse Ising model allows for different types of realizations for the VT- 
symmetry. We easily observe that the macro-reflections can not be implemented 

[P'T,H] ^0, (2.13) 

whereas all the micro-reflections can be realized 

[VT, H] = for k x , K y e iR, [V x / y T, H) = for K x/y G R, K y/x G *R. (2.14) 

Once again these different possibilities raise the question about the unique of the operator 
content in the model. 

Having an explanation for the nature of the eigenvalue spectra, it is left to show that 
one may in addition construct a meaningful metric for this Hamiltonian with well defined 
quantum mechanical observables associated to it. As already indicated, the metric is not 
even expected to be unique so that, unlike as for the Hermitian case, the observables are no 



longer defined by the Hamiltonian alone [ 21 1 . It remains therefore ambiguous what Hamil- 
tonians of the type H (A, k) describe in terms of physical observables. Having constructed 
a metric one may often also compute an isospectral Hermitian counterpart for H(X, n) for 
which the physical observables have the standard meaning. 

One of the main purposes of this manuscript is that of finding the Hermitian coun- 
terparts of the Hamiltonian Ql.l| ) and studying in some detail (at least for small N) how 
many such Hermitian Hamiltonians can be constructed. 

3. Generalities 

3.1 A new metric and an isospectral Hermitian partner from PT-symmetry 

For the sake of self-consistency, we briefly recall the well known procedure [12, 1^, 14, 15, 16] 
of how to construct a meaningful metric and isospectral Hermitian counterpart, h, for a 
non-Hermitian Hamiltonian, H. We assume the Hamiltonian to be diagonalizable and to 
possess a discrete spectrum. Being non-Hermitian the Hamiltonian has non identical left 
|$) and right eigenvectors |\&) with eigenvalue equations 

H\<S> n )=£ n \$ n ) and fft |tf n ) = e n |* n ) for n G N. (3.1) 
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The eigenvectors are in general not orthogonal (<J> n |<3? m ) 7^ 5 nm , but form a biorthonormal 
basis 

J]|* n )($ n | =1. (3.2) 

n 

We assume the existence of a selfadjoint, but not necessarily positive, parity operator V 
whose adjoint action conjugates the Hamiltonian 

if _ VHV with p2 _ L ^3 3^ 

The action of this operator on the eigenvectors 

V\$ n ) = s n \^ n ) withs n = ±l (3.4) 
defines the signature s = (si, S2, ■ ■ ■ , s n ), which serves to introduce the so-called C-operator 1 

C:=^a„|* n )(* n |, (3.5) 

n 

satisfying 

[C,H]=0, [C,VT]=0, C 2 = I. (3.6) 

Next we employ this operator to define a new operator p, which also relates the Hamiltonian 
to its conjugate 

p := VC, Wp = pH. (3.7) 

Depending now on the assumptions made for p, such systems allow for different types of 
conclusions. When p is positive and Hermitian, but not necessarily invertible, the system 
is referred to as quasi-Hermitian [22, 21]. In this case the existence of a definite metric is 
guaranteed and the eigenvalues are real. In turn when p is invertible and Hermitian, but 
not necessarily positive, the system is called pseudo-Hermitian [23, 24, 25 1. For this type 
of scenario the eigenvalues are always real but no definite conclusions can be made with 
regard to the existence of a definite metric. Here we will identify operators p which are 
quasi-Hermitian as well as pseudo-Hermitian. 

Finally we may factorize p into a new operator 2 r\ and use it to construct an isospectral 
Hermitian counterpart for H 

h = r^Hrf 1 = /i f 44> = pHp' 1 with p = rfr]. (3.8) 



In other words assuming the existence of an inverse for p and its factorization in form of 
Q3.Sj ) one can derive a Hermitian counterpart h for H and vice versa. 

1 The is an unfortunate notation and it should be pointed out that the operator is not related to the 
standard charge conjugation operator in quantum field theory. 

2 When 77 is Hermitian, it just corresponds to a Dyson transformation [^6[ employed in the so-called 
Holstein-Primakov method For practical purposes it is useful to have a name for this operator and 

therefore we refer to r\ from now on as the Dyson map. 
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3.2 Expectation values of local observables 

As discussed above, when dealing with non-Hermitian Hamiltonians the standard metric 
is generally indefinite and therefore a new, physically sensible, metric needs to be defined 
by means of the construction described before. This amounts to introducing a new inner 
product ( | ) p which is defined in terms of the standard inner product ( j ) as 

(*|*) p := (*|ptf>, (3.9) 

for arbitrary states, ($| and Assuming that all local operators O in the non-Hermitian 
theory are related to their counterparts o in the Hermitian theory in the same manner as 
the corresponding Hamiltonians 

rjOn' 1 = o, (3.10) 
one finds that a generic matrix element of the operator O has the form, 

($|p0|tf> = ($\ri*on\V) = (<£|e#> , (3.11) 

where \^>) and (<I>| are eigenstates of the non-Hermitian Hamiltonian and its conjugate, 
respectively. The states and {<f>\ are related to the previous two states by = r]\^) and 
= (3>|^, that is, they are eigenstates of the Hermitian Hamiltonian corresponding to the 
same eigenvalues. Equation ( |3.11| ) will be used later on in this paper for the computation 
of various kinds of expectation values. 

3.3 Perturbation theory 

In most cases the above mentioned operators can not be computed exactly and one has 
to resort to a perturbative analysis. Let us recall the main features of such a treatment. 
To start with it is convenient to separate the Hamiltonian into its Hermitian and non- 
Hermitian part as H(X,n) = ho(\) + inhi, where ho and hi are both Hermitian with k 
being a real coupling constant. The latter term may then be treated as the perturbing 



term. For our concrete case (1J.) the individual components are 



N N 

MA) = ~2 J>? + A ^<+i)' and fci = -2E°?> ( 3 - 12 ) 

8=1 1=1 

such that ho(\) corresponds to the Ising spin chain coupled to a magnetic field in the z 
direction and the perturbing term is an imaginary magnetic field in the x-direction. In 
order to determine n, p and h we can now solve either of the two equations in ([T^). Here 
we decide to commence with the latter. Making the further assumption that n is Hermitian 
and of the form rj = e q l 2 this amounts to solving 

fft = e 9 He -i =H+[q,H] + ~[g, [q, H]] + ^[q, [q, [q, H]}] +■■■ (3.13) 

where we have employed the Backer-Campbell-Hausdorff identity. Writing H and in 



terms of ho and hi equation (3.13) becomes 



%nh x + %K[q, hi] + y [q, [q, hi]] +■■■ = [ho, q] + -[q, [h , ?]]+■•■ (3-14) 
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For most non-Hermitian Hamiltonians, such as for our model ( |1.1[) , this equation is very 
difficult to solve for q. When the (^+l)-fold commutator of q with ho, denoted by c{ +l \ho) 
vanishes, closed formulae were found in [28] 

2 ( l\ n 7? 2 

h = h » + E lW cf " )(/lo) ' H = ho ~ £ (^T)I c ? n " 1)(/lo) ' (3 - 15) 

where [x] denotes the integer part of a number x and E n are Euler's numbers, e.g. E\ = 1, 
£/2 = 5, £3 = 61, -E4 = 1385, . . . The coefficients «2n-i were determined by means of a 
recursive equation, which was solved by 

, [(n+l)/2] , v 

m=l v 7 

such that ki = 1/2, k 3 = -1/4, k 5 = 1/2, k 7 = -17/8, . . . 

One may also impose some further structure on q and expand it as 

00 

! = E^"Wi, (3-17) 
k=i 

so that each perturbative contribution q^k-i is a ^-independent matrix. For models of the 
form considered here only odd powers of k appear in the perturbative expansion. This is 
essentially due to the fact that H and W are related to each other by n — > —re. Substituting 
the expansion ( [3.17 ) into the equation ( |3.14j ) one finds a set of equations for q\, q$, q$, . . . 



by equating those terms in ( |3.14|) which are of the same order in perturbation theory in re. 
The first few equations are given by 

[h ,q 1 ]=2ih 1 , (3.18) 
[ho,q 3 } = ~[gi,[<7l,ML (3-19) 

% % % 

[ho,Qs] = gfaii [Q3,hi]] + -[q 3 , [quh]] - ggpfeli [91, [ft, (3-20) 

As we can see easily, they can be solved recursively, namely once q\ is known, one case solve 
for (73 and so on. A closed expression for the commutator [ho, q n ] in terms of commutators 



[g m ,/ii] with m < n was derived in [18|. Perturbation theory has been carried out in the 



past for various non-Hermitian models, e.g. |16|, ^8[ |2^, g(J 31], UI- 

The model at hand is special in the sense that it involves two coupling constants, 
i.e. re and A, such that it allows for an alternative perturbative expansion in terms of the 
latter. Indeed we will demonstrate below that the case A = can be solved exactly and we 
can therefore expand around that solution. Proceeding similarly as for the re-perturbation 
theory we separate the Hamiltonian into its single spin contribution and into the nearest 
neighbour interaction term H (A, re) = Hq(k) + A/iiwith 

N N 

i=l i=i 
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We stress that the counterparts of ( |3.18| )- (|3.20| ) in the well known K-expansion explained 
above differ substantially in the A-expansion. The details will be explained in the main 
part of the manuscript below. Having the option to construct two perturbative series, we 
in principle have in addition the possibility to combine them in a manner that has proved 



to be very successful in the context of high intensity laser physics [ 32 



3.4 Ambiguities in the physical observables 

As mentioned previously, one can argue that the metric p is not unique. In the perturbation 
theory framework, this can be easily seen from the fact that the equations ( 3.18 )-( j3.20|) (and 
any other equations arising at higher orders in perturbation theory) admit many different 
solutions. The non-uniqueness of r\ or, equivalently, the fact that several independent 
Hermitian Hamiltonians h may exist which are all related to the same non-Hermitian 
Hamiltonian by different unitary transformations is well known in the literature. Indeed, 



34, 35, 36] and is currently still 



this fact has been noticed already in the past [21, 33, 
object of debate [37, |3§| ]. 

Assuming now the Dyson map rj in fl3.S| ) to be Hermitian and related to the operators 
V, C and p as defined in (|3.7D we simply obtain 



rj = ri ] rj 2 = p = VC. (3.22) 

Writing 77 = e 9//2 , it is obvious that we can always add to q any matrix b that commutes 
with the full Hamiltonian [H, b] = and with q, [q, b] = 



h 



(3.23) 



and still solve equations ( |3,8j ) . This kind of ambiguity is not very interesting, as it will not 
change h and therefore not lead to new physics. A somewhat less trivial ambiguity was 



pointed out in [30|, which will generate different types of Hermitian counter-parts to H. It 
originates from the fact that we can always add to q%, q%, #5, . . . any matrix commuting with 
ho as we may easily observe in equations (|3.18| )- (|3.20 ). Below we will see that in principle 
for specific examples many such matrices can be found. 

However, by relating rj to the operators C and V as in ( |3.22 ) we are introducing further 
constraints on the form of rj. These constraints follow from the equations (|3.6| ), particularly 



the last two equations there. Using the explicit form (3.22) they can be rewritten as 

VTe q VT = e g , Ve q V = e~ q . (3.24) 
r] 2 V 



by employing the equality C 
that 



= e q V. In order for (3.24) to be satisfied, it is required 
VqV = TqT = -q, (3.25) 



and consequently 



Vq 2 k~iP = Tq 2 k-iT = -q 2 k-l, V k G 



(3.26) 



Below, we will show that these constraints are sufficient in many cases to fix the operator r\ 
and therefore the metric completely. However, it should be noted that these arguments are 
based on the assumption that p acquires the form ( 3.22j ) and furthermore that the parity 
operator is unique, which as we exemplified (2.8) is not always the case. 
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4. The Yang-Lee quantum chain: non perturbative results 

We will now employ the general ideas and definitions introduced in the previous subsection 
for the quantum spin chain Hamiltonian (|1.1|). In particular, we will show how to obtain 
exact solutions for the operators rj, p and h in the two particular situations: i) A or k are 
vanishing and N is generic and ii) A and k are arbitrary and N is taken to be small. 
For large values of N it will be convenient to use the following abbreviation 

N 

s a ia2 ...a p ^2 a k la k+i--- a k+p-v for ai = x,y,z,u; i = 1, . . . ,p < N. (4.1) 

k=l 

We denote here a u = I to allow for non-local, i.e. not nearest neighbour, interactions. In 
this notation the Hamiltonian ( |1 . 1|) reads 



H(X,k) = h (\) + iKhx, with h (\) = --(S? + \S* X ), h x = -~S™. (4.2) 

In what follows it will also be important to use the adjoint action of V, T and VT on the 
generators S^ ia2 ^ ap . It is easy to compute 

VS^^V = (-l) n » +n *S» a2 ... ap , (4.3) 

n-qN T—(_-\\ n vQ N (A A\ 

2 u a\a2...a p A \ ^ a\a^...a v i v**^v 

VTS^ a VT = (-1)^5^...^, (4.5) 

where the numbers of indices at equal to x, y, respectively. These identities 

follow directly from the definitions (|2.4|) and (|2.5|). 



4.1 Limiting cases: A = or k = 

Let us start by considering the special case A = for which 

ho(0) = -±S? and h x = ~S%. (4.6) 

Although the Hamiltonian is extremely simple, it is still non-Hermitian, and thus serves 

as a benchmark to illustrate the above mentioned notions. For example, a matrix r) that 
relates H (0, k) to its Hermitian counterpart h(0, k) is easily found to be 

77 = e <i/ 2 = g-|arctanh(/t)S^_ Uj\ 

Its adjoint action on S% and is simply 

nS?^ 1 = - 7 =L= (ikS? + S» ), r^V 1 = 7f=f " iKS *)> ( 4 - 8 ) 



which when we evaluate (3.8) yields the Hermitian counterpart to Hq(k) in ( p.21| ) 

h(0,K) = ~Vl-K 2 S?. (4.9) 
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This Hamiltonian describes a spin chain for which no mutual interaction between spins 
along the chain occurs. An external magnetic field is applied at each site of the chain, 
whose intensity is governed by the value of k and is the same at every site. The constraint 
— 1 < k < 1 ensures the Hamiltonian /i(0, k) and r\ to be Hermitian. Given the simplicity of 
/i(0, n) we can easily find its full set of eigenstates and eigenvalues, hence those of H(0, n). 
The operator is a diagonal matrix with entries 



di&g(N,N-2,...,-N + 2,-N). (4.10) 



The entries in the diagonal (eigenvalues) are N — 2p with p = 0, . . . , N. They are not neces- 
sarily in decreasing order and, except for N and — N, all other eigenvalues are degenerate. 
For example, the eigenvalues N — 2 and 2 — N are always N times degenerate. This means 
that there is a single ground state with minimum energy, 

N 



1 V I 

E 9 (k) = --V1-^, (4.11) 
and the corresponding eigenstate is simply 




associated to a configuration with all spins "up", hence aligned with the magnetic field 
that is being applied at each site of the chain. 

The situation when k = and A is arbitrary corresponds to the Hermitian Hamiltonian 
given by /to (A), that is the Ising spin chain with a magnetic field in the z-direction. In 
this case, 77 = I, which is automatically ensured when using perturbation theory. The 
eigenstates and eigenvalues of this Hamiltonian have been studied in the literature by 



using the Bethe ansatz approach, see e.g. [|39|, 40]. In particular, the ground state can not 



be written in such as simple form as ( 4.12j ), as it will depend on the value of A. One does 



know however, that, for finite N, it will interpolate between the A = case, in which the 
ground state is ( 4,12; ) and the A — » 00 case, in which the ground state will correspond to 



alternating up-down spins. 

4.1.1 Uniqueness of the Dyson operator 



In light of the discussion in section [L4J it is also interesting to investigate the uniqueness of 
J4.7p . Indeed, we will now show that ( |4.7| ) is the only solution to ( |3. 8| ) which is consistent 
with ( |3.240 for the Hamiltonian H(0,k). This can be proven in two steps: firstly we will 



characterize the subset of matrices and linear combinations thereof that satisfy ( 3.25| ) and 



secondly, we will show that none of these matrices can be in the kernel of ho(0). Let us 
define the matrices, which provide a basis for the set of 2 N 2 Ar -Hermitian matrices, 

M ai ...a N = <?T <8> • • • <8> <7^f, with o« = x, y, z, or u Vi = l,...,N. (4.13) 

Recall the definition af = Ij. Let us consider an arbitrary linear combination of the 
matrices ( 4.13j ), The action of parity and time reversal on such a linear combination is 
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analogous to ( t4,3| ) and (|4.4|). From this it follows that, in order for any linear combination 
of matrices M i ) ...,o«, to transform as q does in equations ( 3.25[ ) it must be such that for 



all matrices in the linear combination n y is odd and n x is even (n x and n y as defined after 
equation (|4.5| )), 

We will now argue that no matrix in the kernel of ho(0) is of this form. There are 
various ways of having a vanishing commutator [ho(0), B] = 0. The most obvious solution 
is for B to be a diagonal matrix, as ho(0) is itself diagonal. In terms of the matrices ( 4.13j ), 



this means selecting out those that are tensor products of a z and I only. There are overall 
2^ such matrices and obviously none of them has n y odd. This would be sufficient to 
conclude that the solution ( |4.7| ) is unique if only the kernel of ho(0) had dimension 2^. 
This is not so because ho(0) has degenerate eigenvalues. 

Any additional matrices in the kernel will be some linear combination of matrices 



(4.13) involving at least one index x or y. Employing the commutation relations (|1.3j), it is 
easy to see that there are basically two kinds of additional matrices that are in the kernel 
of ho(0): firstly, the matrices M xyu ___ u — M yxu __ M and generalizations thereof , which are 
antisymmetric under the exchange of indices x <-> y and violate the condition n x even and 
secondly, the matrices M xxu ___ u + M yyu _ u and generalizations thereof, which violate the 
condition n y odd and are symmetric under the exchange of indices x *-* y. Generalizations 
of these matrices are those obtained by replacing any number of indices u by z and/or 
permuting indices, as well as other matrices of similar characteristics, such as M xxxxu _ u + 
M yyyyu ... u + M xxyyu _ u + M yyxxu ^ M and so on. Since this is more an argument than a proof, 
we would like to support it with two examples. For N = 2 

MO) =diag(-l, 0,0,1), (4.14) 

and the kernel has dimensions 6, as one eigenvalue is twice degenerate. It is generated by 
the matrices 

M xy -M yx , M xx + M yy , M zz , M zu , M uz and M uu = I. (4.15) 

For N = 3 we have that /io(0) has four different eigenvalues, two of which are three times 
degenerate, 

, ( 3 1 1 1 1 1 1 3\ 

fto(0) = diag , (4.16) 



2 2 22 2 2 2 2, 
The dimension of the kernel then becomes 20. Its generators are the matrices 

M xyu M yxU} M xuy M yux , M uxy M uyx , 

M xxu -\- M yyu , M xux + M yuy , M uxx + M uyy , 

M xyz M yxz , M xzy M yzx , M zxy M zyx , 

M xx z + M yyz , M xzx + M yzy , M zxx + M zyy , 

M zzz , M zzu , M zuz , M uzz , M zuu , M uzu , M uuz , M uuu , (4.17) 

As shown before, these examples confirm once more that no element in the kernel of ho(0) 
can fulfill the conditions ( 3.25| ) and therefore could not be added to q, whilst fulfilling such 



conditions. Thus no matrices in the kernel of /io(0) satisfy the conditions ( 3.25| ) and the 
solution Q4.7D is unique if the operator r\ = e q l 2 is to be Hermitian. 
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4.2 The N = 2 case: two sites 

We have already identified the PT-symmetry for the Hamiltonian (|1,1| ) with V given as 
specified in fl2,3| ) satisfying (|3.3| ), Let us now take the length of the spin chain to be N = 2 
and compute the quantities as outlined in the previous section. 

For two sites we may chose without loss of generality the boundary conditions to be 
periodic c^+i = a i as an y other choice may be achieved simply by a re-definition of A. In 



this case the Hamiltonian (1.1) acquires the simple form of a non-Hermitian 4 x 4-matrix. 
In order to make notations clear, we will write this matrix here in the various notations 
introduced so far, 



^ [a z <g> I + I ® a z + 2\a x <g> a x + w (I ® a x + a x <g> I)] , 



--[at + a z 2 + 2Xafa x + m {a x 2 + af] 



/ -i IK IK 

I 1 2 2 

A 

A 

\ IK %K 

\ A T T 



A\ 



2 

in 
2 



Hi. 

2 
in 
2 



(4.18) 



where the first line shows the most explicit way of writing the Hamiltonian, the second line 
shows a simplified version, were the tensor products are omitted and absorbed into the as 
as specified after ( |1 . 1[ ) . The last line uses the notation introduced in (p~l|). 



■ 1 1 1 1 1 1 








- u 




P7 







Figure 1: Domains of broken and unbroken 'PT-symmetry 

At first we shall be concerned with the spectral properties of this Hamiltonian. The two 
subdomains U-pq- and Ub-pj- , as introduced in (2.7), have already been identified numerically 
in H for spin chain lengths up to N = 19, that is for matrices up to the remarkable size 
of 524288 x 524288. For N = 2 the eigenvalues for ( 4.18j ) are easily computed analytically 
as the characteristic polynomial factorizes into a third and first order polynomial. The 
discriminant A of the third order polynomial is computed by 

A = r 2 - q 3 with q = ^ (-3k 2 + 4A 2 + 3) , r = A (18k 2 + 8A 2 + 9) . (4.19) 
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The eigenvalues are guaranteed to be real when the discriminant is smaller or equal to zero, 
such that U-pr is defined as 

U VT = { A, k : A = k 6 + 8AV - 3k 4 + 16A 4 k 2 + 20AV + 3k 2 - A 2 - 1 < 0} . (4.20) 

The regions U-pq- and U^-pq- are depicted in Figure 1, from which we note that in order to 
have a real eigenvalue spectrum k is restricted to take values between and 1, whereas A 
is left unbounded A G [0, oo). 
The four real eigenvalues are then computed to 



ei = A, e 2 = 2q\ cos (§) - |, e 3 ,4 = 2q\ cos (f + tt t ~ i 



(4.21) 



where the additional abbreviation 6 = arccos (r/q 3 ^ 2 ) has been introduced. We depict 
these eigenvalues in Figure 2, 





Figure 2: Avoided level crossing: eigenvalues as functions of A (n) for fixed n (A). 

where we observe the typical avoided level crossing behaviour of the eigenvalues as a func- 



tion of the parameters 41], i.e. the eigenvalues £3 and £4 only meet in the exceptional 
point when they simultaneously become complex. 

For the computations of physical observables, which we will carry out below, it is 
important to identify the lowest eigenvalue, which turns out to be always £4. 

Next we compute the right eigenvectors of H (A, k) to 



|$a) = (0,-1,-1,0), |* n ) = (f n ,-a n ,-a n ,P n ), n = 2,3,4, 



(4.22) 



with a n = ik (A — £ n + 1), n = k 2 + 2A 2 + 2A£ n and j r 



2el + 2A - 2A£ n + 2£ r , 



We verify that left and right eigenvectors are related via a conjugation |^ n ) = ($ n | and 
compute the signature as defined in ([O]) to s = (+,—,+, — ) for the parity operator (fO|). 
Normalizing the vectors in Q4.22[ ) by dividing with Ni = \/2, N n = (2a 2 + /3 2 + 7 2 ) 1 / 2 for 
n = 2, 3, 4 we compute the C-operator according to (p.5[) to 



C 



( c 5 


-c 3 


-c 3 






-c 3 


-Ci-1 




c 2 




-c 3 


-Ci 


-Ci - 


1 c 2 




V c A 


c 2 


c 2 


2(Ci + 1) 


-c 5 ) 



(4.23) 
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where the matrix entries are 



r, — fti"^ 



A'.V 



<4 



1 

2' 



+ 



0373 _ Pili 



c 2 
c 5 



7l 



At 



+ 



aT 



Nf 



"272 1 "373 



0474 



We may now verify that C indeed satisfied the properties (]3_. 

G2 — O2O5 — O3U4, O3 — O5U3 — O2O4 — ZU1O3 

1 = 2Cf + C| + Cf , = C| + Cf + 2Ci (Ci + 1). 
Next we compute the metric operator in the form p = 



(4.24) 

upon the use of the identities 

4 — O2O3 — C1C4, . . 

(4.25) 

VC simply from (|]| and (|4.23| ) 



to 



/C 5 -C 3 
C 3 1 + d 
C3 Ci 

\ C4 C2 



— C3 C4 

Ci — C2 

1 + Ci -c 2 

C 2 2(1 + Ci)- 



\ 



(4.26) 



c 5 y 

Since ioi, /9j,7j S 1 it follows that Ci, id, id, C4, C5 £ R and therefore we conclude 
immediately that p is Hermitian. To see whether p is also positive, as it ought to be, we 
compute its eigenvalues 



y 1= y 2 = l and y 3/4 = 1 + 2Ci ± 2y/C 1 (l + C 1 ). 

Since Ci > all eigenvalues of p are obviously guaranteed to be positive. 
Next we determine the corresponding eigenstates to 



(4.27) 



Fl 



(0,-1,1,0), |r 2 ) = (C 4 , 0,0,1 



C s ), 
- l))/2 



1^3/4} = (73/4> "3/4, "3/4, ^3/4) (4-28) 



with a 3/4 = y 3/4 (C 3 C A + C 2 (-4d + C 5 - l))/2 - C3C4, £3/4 
(Cf + d (4Ci - C 5 + 3)) y 3/4 and 7 3/4 = C1C4 - C 2 C 3 + (C 2 C 3 + C 1 C 4 )y: 
the matrix U = {ri, r 2 , r%, r 4 }, whose column vectors are the eigenvectors of p, we may 



C\ — C1C5 + 
3/4. Defining now 



take the square root of p, such that r\ 



= UD l l 2 U \ where D = diag(yi,y 2i , y 3) , y 4 ). 

to an spin chain (with 



The isospectral Hermitian counterpart of H results from (3. 
just two sites) in a magnetic field 

h{X,n) = r/tfrT 1 = UD 1/2 U- 1 HUD- 1/2 U- 1 

P 2 zz (X,k)S 2 zz + p 2 z (X,n)S 2 



P 2 xx (X,k)SI x + p 2 (X,k)S 2 



yy^ 
,2 



(4.29) 
(4.30) 



L HUD- 

\, ^Syy 

It is clear that the coefficients p xx , p 2 y , p 2 z , p z can be computed explicitly, but the 
expressions are rather lengthy and we will therefore not present them here. They are all real 
functions of A and k. Their explicit form can be found in appendix |A| in terms of functions of 
A and k fl5.10j ) which will be introduced in section ||, in the context of perturbation theory. 
In the next section we wish to compare this exact result with a perturbative expansion. 
Let us therefore report two numerical examples for some isospectral Hermitian counterpart 
of H(X, k) 

( -0.829536 -0.0606492 \ 

-0.0341687 -0.1341687 
-0.1341687 -0.0341687 
\ -0.0606492 0.897873 / 



/i(0.1,0.5) 



(4.31) 
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/i(0.9,0.1) 



(4.32) 



and 

/ -0.985439 -0.890532 \ 

-0.0094167 -0.909417 
-0.909417 -0.0094167 
\ -0.890532 1.00427 / 

Notice that h 2 3 = ^32 = ^22 — A = /133 — A. 

We have carried out a similar analysis for the chain with three sites explicitly, albeit 
the resulting formulae are rather cumbersome to present. In any case for longer chains 
one has to resort to more sophisticated and less transparent techniques as for instance the 
Bethe ansatz. Alternatively, we may employ perturbation theory. 



5. The Yang-Lee quantum chain: perturbative results 

In this section we want to address the problem of obtaining the matrices rj, p and h from 



a perturbative analysis as described in section 3.2. We will study the N = 2, 3 and 4 cases 
in detail and draw some conclusions concerning the analytic expressions of 77, p and h for 
generic N. 

5.1 The N = 2 case: perturbation theory in k 

Despite the fact that H(X, k) is just the 4x4-matrix ( 4.18] ), it is actually not easy to find the 



matrix q in ( 3.13 ) exactly. As discussed in section |3.4| , it is clear that the equations ( 3.18j ) 
to ( [3. 20 ) as well as the equations that would be obtained for higher orders in perturbation 
theory, admit many solutions. Any solution g2fc-l can be modified by adding a matrix that 
commutes with /10(A). However, not all solutions obtained in this manner would be valid 
solutions if the equations (3.24) are to hold. For the particular case N = 2, we are about 
to show that these constraints actually select out a unique Hermitian counterpart to the 
Hamiltonian H(X,k). We will start by finding the most general matrix qi(\) which solves 
the identity (|3.18| ). It is quite clear that given one solution qi(X), any matrix of the form 
qi(X) + B(X) with [ho(X), B(X)] = will also be a solution, so we may start by finding 
all such matrices. In this simple case, there are four basic independent solutions to the 
equation [ho(X), B(X)] = 



i?i = I, B2 = S, 



q2 1 ql 

xx ' yy 



and B4 = Si — AS 1 .; 



mi ■ 



(5.1) 



Since Hq{X) is a 4 x 4-diagonalizable matrix, with non-degenerate eigenvalues, there can 
be at most four independent matrices that commute with it, namely those shown above or 
combinations thereof. On the other hand, it is clear that any polynomial function of the 
Hamiltonian Kq{X) would also commute with (A). As the four matrices in ( p^ ) constitute 
a basis, we expect to be able to express any power of ho as linear combinations of them. 
Indeed, we find 



MA) 2n 

MA) 2n+1 



(1 + A 2 )" 1 . A 2 ™ 
(B l + -B 2 ) + - 



(B 



1 



2\n 



(i+A 2 rMA) + 



A(l + A^) 



2 

^2n+l 



Bo 



-Btk 



(5.2) 
(5.3) 
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for n 6 No. Therefore, the most general solution to the first order equation ( 3.18 ) for the 
present model is 

4 

qi (\) = -Si - X(S 2 yz + S%) + Y, fiW B i, (5-4) 

i=i 

where the fi(\), i = 1, 2, 3, 4 are arbitrary functions of A. 

Before we proceed to determine (/3(A) by solving ( 3.1S| ) let us comment on the am- 
biguities and answer the question of whether all solutions ( |5.4|) are compatible with the 



equations (|3.24[) . Specializing equations flO| ) and (|4.4[) for the matrices in (|5.4|) we find 

VXV = TXT =-X, for X = Sj, S 2 Z , S 2 zy (5.5) 

whereas 

VBiV = TBiT = B h for i = 1,2,3,4. (5.6) 

These equations imply that the equalities fl3.24|) can only be satisfied if the functions 
fi(X) = for i = 1,2,3,4. Thus we have selected out a unique solution for (71(A), namely 

qi (X) = -S 2 -X(S 2 yz + S 2 zy ). (5.7) 

More generally, the conditions ( |3,24| ) together with the properties ( |4.3| ) and ( |4,4| ) imply 
that 

• any solutions (72/fc— 1 must be linear combinations of matrices ( |4.1| ) with odd, 

• any solutions (?2fc— 1 must be linear combinations of matrices ( |4.1| ) with n y + n x odd, 



• or, combining the two conditions above, any solutions (fejfc-i must be linear combina- 
tions of matrices ( |4.1| ) with n y odd and even, 

as anticipated in subsection |4.1.1| . These conditions then automatically guarantee the 
validity of the "PT-properties ( |3.26| ) for the ?2fc— l- F° r N = 2, this singles out the matrices 
S 2 and S yz = S zy in (5.7), so that, even before attempting to solve ( |3.18 ) we would already 
know that it can only be a linear combination of those two matrices. As indicated above, 
these constraints apply for all other q^fc— l(A), with k > 1 so that we can safely claim that, 
at all orders in perturbation theory, the matrices (?2fc-i(A) must be linear combinations of 
the form, 

g 2fc -i(A) = a 2k ^(X)S 2 + b 2k _ 1 (X)(S 2 yz + S 2 zy ), (5.8) 

where 02^-1 (A), &2fc-i(A) are real functions of A. In other words, all the terms in the 
perturbative expansion of q are linear combinations of the same two matrices. Hence, we 
can write 



e <? = e a(X,K)S^+(3(\,K)(S^ z +S 2 zy ) ^ ( 5 _g) 



which, after computing the exponential becomes 

/ p(A,k) 2 +<e(A,k) 2 cosh[27(A,/c)] ie(A,/t) sinh[27(A,«)] ie(X,n) sinh[27(A,K)] S(X,k) sinh 2 [7(A,ftQ] \ 

27(A,k) 2 27(A,k) 27(A.k) 7(A,k) 2 

" (A ' K 27ff) 7(A ' K)1 cosh 2 7 (A, K ) sinh 2 7 (A,K) _ i P (\^pj{X, K )] 

ie(X,K)smh[2sy(X,K)] sinh 2 7 (A, k) COsh 2 7 (A, k) _!P&|gM 

<5(A,k) sinh 2 7(A,k) ip(A,K) sinh[27(A,K)] ip(X,n) sinh[27(A,«;)] <e(A,k) 2 +p(A,k) 2 cosh[27(A,K)] 

\ 7(A,ft) 2 27(A,k) 27(A,k) 27(A,k) 2 
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where 

oo oo 

a(X,K) = 5> 2fe+1 a 2fc+ i(A), P(X,k) = £ K 2k+1 b 2k+1 (X), (5.10) 

k=0 k=0 

and 

7 (A, k) = V«(A, k) 2 + 4/?(A, k) 2 , S(X, k) = a(A, k) 2 - 4/3(A, k) 2 . (5.11) 

e(A, «) = a(A, k) + 2/3(A, «), p(A, k) = a(A, «) - 2/?(A, k). (5.12) 

Notice that, for a(X, k) and (3(X,k) real, the matrix above is explicitly Hermitian, as it 
should be. Once the coefficients ck(A,k) and /3(A,k) have been obtained, the Hermitian 
Hamiltonian ([T^) can be easily computed. The difficulty here is however that general 
formulae for the coefficients a2fc+i(A) and b 2 k+i (A) are very difficult to obtain. Nonetheless, 
perturbation theory allows us to compute these coefficients up to very high orders in powers 
of k. In order to solve for such high orders, we have resorted to the use of the algebraic 
manipulation software Mathematica. It allows us to find the entries of the matrix ( |5.9| ) as 
perturbative series in k and to fix the coefficients 02^+1 (A) and b 2 k+i(X) by matching the 
entries of W(X,k) and r] 2 H(X, k)t]~~ 2 , order by order in perturbation theory, as expected 
from ( |3.8| ). For numerical computations and sufficiently small values of k this gives results 
which are very close to the exact values. In tables 1 and 2 we present the coefficients 
a2fc+i(A) and b 2k+ i(X) up to k = 7. 





-A u 


-A 


-A 


-A" 


-A 


-A 


-A 


-A 


ai(X) 


1 























as (A) 


1 

3 


2 4 
3 




















015(A) 


1 

5 


244 
15 


2* 

5 

















07(A) 


1 

7 


1152 
35 


35104 
105 


2 1 ' 1 
7 














09(A) 


1 

9 


17432 
315 


43408 
35 


1890368 
315 


old 

IT 











011 (A) 


1 

11 


289616 
3465 


797296 
231 


38228224 
1155 


355526144 
3465 


2™ 
11 








013(A) 


1 

13 


353372 
3003 


72293440 
9009 


655729408 
5005 


2275245568 
3003 


15442769920 
9009 


2 M 
13 





015(A) 


1 

15 


7100416 
45045 


67453952 
4095 


896579072 
2145 


58903814656 
15015 


717363822592 
45045 


1273503367168 
45045 


15 



Table 1: The coefficients a 2 fc+i(A) for ft < 8. 
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-A 


-A 


-A 


-A 


-A 


-A 


-A ib 
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63(A) 


4 
3 


2 4 
3 




















65(A) 


23 
15 


2' 
5 


t. 

















07(A) 


176 
105 


7544 
105 


2 lJ (3) 
7 


2 12 

7 














09(A) 


563 
315 


49136 
315 


212816 
105 


9 


2 m 
9 











611(A) 


6508 
3465 


335576 
1155 


7827328 
1155 


164005504 
3465 


2 la (5) 
11 


2 20 

11 








613(A) 


88069 
45045 


4400960 
9009 


39578944 
2145 


1947324416 
9009 


9037578752 
9009 


2^(3) 
13 


2 24 

13 





615(A) 


91072 
45045 


34381136 
45045 


178162048 
4095 


1068366848 
1365 


37570428928 
6435 


903387164672 
45045 


2 AB (7) 
15 


2 28 

15 



Table 2: The coefficients o 2j t + i(A) for ft < 8. 



These tables should be understood as follows: in order to obtain the corresponding coef- 
ficient the numbers in a given row are to be multiplied by the power of A (with a minus 
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sign added) at the top of the same column and added up. For example: 

1 244A 2 2 8 A 4 

" 5 < A > = -5--nr-— (5 - 13) 

The only case for which it is easy to conjecture the expressions of a2fc+i(A), &2fc+i(A) for 
generic values of k corresponds to A = 0. Then a2fc+i(0) = —l/(2k + 1) and b 2 k+i(0) = 0, 
which gives the already known result a(0, k) = — arctanh(fv) and f3(0, k) = 0, see section [P| . 
Having found rj, it is straightforward using ( |3.8[ ) to determine the Hermitian counterpart 
of H(X, k). In general, we find 



(h 11 {X,K) h 14 (X,K)\ 

h 22 (X,K) h 22 (X,K)-X 
h 22 (X, K )-X h 22 (X,K) 

yhi4(X,K) /144(A) ft) J 



_ h 22 (X, n)-\ + hu(X, k) 2 h 22 (X, k) - X - h M (X, k) 2 

— ^ b xx + ^ &yy 

| /j-n(A, k) + fi-44(A, k) - 2h 22 (X, k) q2 i hn(X,K) - /t 44 (A,/t) q2 

| /tn (A, ft) + /t44(A, k) + 2h 22 (X, k) ^ ^ 

which is the same kind of structure found in ( 4.3C| ). The functions /in (A, k), h 22 (X,K), 
/ii4(A, k) and /144(A) k) are real functions of the coupling constants which can be evaluated 
very accurately for fixed values of A and k by using the perturbative results above. In 
fact, the remaining entries of the matrix are not explicitly zero as functions of a(A, n) 
and /3(X,k). They are complicated functions of the latter which when carrying out the 
perturbation theory result to be zero up order k 15 . This is consistent with the exact 
results obtained before. The explicit expressions of the entries of h(X, k) in terms of the 
functions ( 5.11| ) and ( |5.12 ) can be found in appendix K[ Here, we will just present their 
expression as a series expansion in k up to order k 4 (for higher orders, expression become 
too cumbersome), 

/in (A, .) = -l + Q + Ap + Q + A + ^ + 4A 3 ) k 4 + 0(k 6 ), (5.15) 
MA, «) = -Ak 2 - A (1 + 4A 2 ) k 4 + 0(n e ), (5.16) 
h M {X, k) = 1 + + X^Jk 2 + ^-i + A - ^- + 4A 3 ^) + 0(k 6 ), (5.17) 

/ii 4 (A, k) = -A + Ak 2 + ( y + 4A 3 ^) k 4 + 0(k 6 ). (5.18) 

From this expansions we can deduce some interesting features which also extend to higher 
orders in perturbation theory 

/ill (— A, k) = —/144(A) k), h 22 (-X, k) = —/122(A) k), hu(—X, k) = -/ii 4 (A, re). (5.19) 

Furthermore, we note that the Hermitian Hamiltonian h(X, re) is an even function of re, 
so that the series expansion of its components involves only even powers of the coupling. 
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Finally, as it should be, the Hamiltonian h(X, k) is also PT-symmetric, which follows from 
the fact that all matrices involved (S%. x , S yy , S% z and S%) are invariant under the adjoint 
action of the operator VT. These are in fact the only matrices that are both PT-symmetric 
and real. In fact we could have known a priori before carrying any computations that h(X, n) 
has to be some linear combination of S% x , S yy , S^ z and S%. Notice that the reality of h(X, k) 

l|) involved must have n y even, as defined 



can be expressed by saying that any matrices 
in the paragraph after equation (jl.5|). 

In order to compare with the results obtained in section [| we give below the numerical 
values of the entries of the Hermitian Hamiltonian /i(A, k) for fixed values of the couplings 



/i(0.1,0.5) 



and 



/i(0.9,0.1) 



/ -0.829534 



\ -0.0606716 



-0.985439 



-0.890532 



V 





-0.0341688 
-0.134169 






-0.00941674 
-0.909417 






-0.134169 
-0.0341688 




-0.0606716 \ 



0.897872 / 



(5.20) 





-0.909417 
-0.00941674 




-0.890532 \ 



1.00427 / 



(5.21) 



We underlined the digits which differ from the exact values computed in ( [4.31 ) and ( 4.32j ) 
and note that the perturbative expressions for h(0. 1,0.5) and h(0.9, 0.1) agree extremely 
well with them, especially for smaller values of k, as is expected. 

In order to see how fast this precision is reached in the perturbation theory we report 
in table 3 the relative error for the entry h\\ order by order up to 15 



X,k\0(k) 


2 


4 


6 


8 


10 


12 


14 


0.9,0.1 


5.7 10" 4 


4.6 10" 5 


4.7 10~ 6 


5.3 10~ 7 


6.4 10~ 8 


8.2 10~ 9 


1.1 10~ 9 


0.1,0.5 


2.5 10~ 2 


6.3 10~ 3 


2.1 10~ 3 


7.5 10~ 4 


2.9 10~ 4 


1.6 10~ 4 


4.7 10~ 5 



Table 3: Relative error = — (perturbative value - exact value) / exact value — for hn order by 
order. 

We observe that the convergence is fairly fast, which allows to extract useful information 
from the perturbation theory even at low order. We shall not be concerned here with more 
rigorous mathematical arguments regarding the summability and convergence in general. 



5.2 The N = 2 case: perturbation theory in A 

In the previous section we have employed the standard version of perturbation theory 
when dealing with non-Hermitian Hamiltonians of the type ( |1.1| ), that is decomposing the 
Hamiltonian into a Hermitian and a non-Hermitian part as in fl3.12|) and then treating 
the non-Hermitian part as the perturbation. Since the Hamiltonian (|1.1| ) depends on two 
independent coupling constants, k and A, it is also natural, albeit less standard, to consider 
perturbation theory in the coupling constant A rather than in k. In other words we expand 
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around the exact solution for A = provided in section 4.1 and treat the nearest neighbour 
interaction term as perturbation. As announced already in section 3.2., we decompose 
H(X, k) into 

H(X,K) = H { K ) + Xh 1 , where H (k) = ~(S? + wSj?), h x = -\s% c . (5.22) 
We wish now once again to solve the equations ([T^) for the Dyson map r], that is 

H\\, k) = e w H(X, K)e~ w , (5.23) 
where we have assumed that n admits the exponential form 

rj = e w/2 with w = ^A a u> a ( K ). (5.24) 

a=0 

At order A equation ( 5.23j ) becomes simply 

Hl( K ) = e wo{K) H (K)e- wo(K) . (5.25) 



The solution to this equation for all ./V was found in subsection 4.1 and corresponds to the 



Dyson map identified in equation (|4.7| ). For N = 2 this means that 

Wq(k) = — arctanh(«:)S'y . (5.26) 

Employing the once again the Backer-Campbell-Hausdorff identity to select 0(A) terms in 
Q5.23 ) we find the condition 



00 k _. 

h 1 = e w °Wh 1 e- w °W+J2Y, E uKW-KW,- ,[w ah (K),H (K)]---]] (5.27) 



k=l i=l Oj=l, 



Notice that, because of the presence of the zeroth order term wq(k), the equation ( 5.27| ) 



involves a sum of infinitely many contributions, as would equations corresponding to higher 
orders in perturbation theory. Because of this, it would in general be difficult to solve ( 5.23|) 



using perturbation theory in A. However, for N = 2 we can solve up to high orders in A by 
exploiting the fact that r\ must have the structure identified in the previous section. This 
means that rj is a matrix of the form ([T^) with 

00 00 
a(X,K) = ^2x a y a ( K ), P(X,k) = ^A a z a ( K ). (5.28) 

a=0 a=0 



It is then possible to find the real functions y a {^) and £ (k) which solve equation ( 5.23] ) 



order by order in A by employing Mathematica, as explained in the previous subsection. 
In this way, we have obtained the functions y a {^) and z a {n) above up to order A 15 . Here 
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we will just report the first five orders, 

Vo («) 
zi(«) 



— arctanh(K), 




(5.29) 


yo(«) 
i-k 2 ' 




(5.30) 


2(k + 2k 3 + (1 - k 2 


) VOW) 


(5.31) 


(1-k 2 ) 3 




2(k + 2k 3 + (1 - k 2 


- 2k 4 ) y ( K )) 


(5.32) 


(1 - k 2 ) 


4 ' 


2 (k (3 - 5k 2 - 32k 4 - 


- 8k 6 ) + (3 - 6k 2 - 5k 4 + 8k 6 ) y o («0) 


(5.33) 




(1 - k 2 ) 6 


2 (k (3 - 5k 2 - 36k 4 - 


- 16k 6 ) + (3 - 6k 2 - 9k 4 + 28k 6 + 8k 8 ) y (/e)) 


(5.34) 




(1-k 2 ) 7 



3/2 («) 
«3(«) 

y4(«) 

«5(«) 

and y2a+i(^) = ^2a(^) = for all a = 0, 1, . . . From these formulae, it is possible to find 
an expression for the Hermitian Hamiltonian h(X, k) as a perturbative series in A. As it 
should be, one finds the same structure ( 5.14j ) with 

, , , r k 2 A 6(-2 + k 2 + 2^T^)a 2 4k 4 A 3 

/ill (A, Kj = — V 1 — K z + 



l~K 2 (1 _ K 2)f (1- K 2 ) 4 

2 (40 - 44k 2 - 57k 4 + 28k 6 + 8Vl - k 2 (-5 + 3k 2 + 8k 4 )) A 4 



+ 0(A 5 ), (5.35) 



(1-k 2 )^ 
k 2 A 4k 4 A 3 



MA, k) = 2 - -j + 0(A & ), (5.36) 



1-K 2 (1 _ K 2)4 
K 2 A 



6 (-2 + k 2 + 2 V / T r ^ 2 ) A 2 4k 4 A 3 



/i 4 4(A,k) = a/1 - k 2 + H 5 1", o 1 

l-« 2 (1-K 2 )2 (1-K 2 ) 4 

2 f40 - 44k 2 - 57k 4 + 28k 6 + 8\/l - k 2 (-5 + 3k 2 + 8k 4 )) A 4 
+— n ^ + 0(A 5 ), (5.37) 

(1-K 2 )" 



-4 + 4k 2 + 3V1-k 2 ) A 

hu(X,K) 



:i-K 2 )t 



4 (8 - 10k 2 - 2k 4 + 4k 6 + Vl - t 2 (2 + k 2 ) (-4 + 5k 2 )) A 



3 

5 



(1-K 2 )5 



+ O(A ). (5.38) 



Notice that the same symmetries ( |5.19 ) are also found here. We also see once again that 
h(X, k) is an even function of k, as only even powers are involved. Computing again 
numerical values for /i(0. 1,0.5) and /i(0.9,0.1) we find almost perfect agreement with the 
exact results. There is extremely good agreement both with the exact results ( [1.31 ) and 
(432) and with the result from perturbation theory in k flpOD and |[2g ). In order to 
see how fast this precision is reached in the perturbation theory we report in table 4 the 
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relative error for the entry hu order by order up to 15, omitting the odd orders despite the 
fact that they occur in the A-perturbation theory 



X,k\0(X) 


2 


4 


6 


8 


10 


12 


14 


0.9,0.1 


3.4 10~ 3 


2.3 10~ 5 


1.9 10~ 6 


1.8 10" 7 


1.9 10~ 8 


2.0 10~ 9 


2.2 10~ 10 


0.1,0.5 


1.1 10~ 3 


6.3 10~ 5 


4.9 10~ 6 


3.6 10" 7 


3.1 10~ 8 


2.8 10~ 9 


2.6 10~ 10 



Table 4: Relative error 
order. 



-(perturbative value - exact value) / exact value — for hu order by 



We note that the perturbation theory converges extremely fast, even for large values of A, 
for which one would not expect such a behaviour. This can be explained as follows: In the 
domain of unbroken PT-symmetry U-pj- the allowed values for k become very small as A 
increases. As we note from the expressions fl5.35|) -( |5.38] ) the order of k increases with the 
order of A term by term. 

5.3 The N = 3 case 

We will now carry out an analogous perturbative study in k for the three sites case. We 
keep the choice of periodic boundary condition, even though for sites more than two this 
means some loss of generality. Proceeding as before, we will try to obtain the matrix q 
perturbatively, by solving the consistency conditions ( |3.18|) - (|3.20|) . Now we have to solve 
the problem for 8 x 8-matrices. We commence by computing the kernel of /io 

B\ =l, = — XSy yz , 53 = XSyy ~ ^ ~ A^ ^ Syy j, ~ S^g-g , B ± = ~ ^yj., 



B, 



Szzzi -^6 — S xyz — Sy XZ , B 7 — XS XX + S z — —2ho(X), 



q3 , c3 i \ c3 

u xx ' yy ' vvz? 



ill)'-' 



in addition to this eight matrices, there are another four, due to the fact that two of the 
eigenvalues of Hq{X) are degenerate. Hence the dimension of the kernel is 12, 



o»o» + oZ<tZ + o%o5, (5.39) 



B g — S z - X(Syy + S 3 ZZ - a\a\ - a\a\ - of erf), B w — u 2 u 3 t« 2 «3 t u 2"3' 
B u = S zz + Xb xxz - X(a 1 + a 3 + a 1 a 2 a 3 + a\a 2 a\), B 12 = a 3 - a l a 2 a 3 - afa^a 3 , 



with [Bi, ho(X)] = for i = 1, . . . , 12. Similarly as in the case N = 2 we find that all of 
these matrices are parity invariant 



VBiV = B h V 



1, 



(5.40) 



which from equations ( |3.26 ) means that no linear combination of the matrices Bi can be 
added to q2k-i that would be compatible with the constraints ( |3.24j ). Therefore, with such 
constraints, there is a unique solution to ( 3.18j ) which has the form, 



gi (A) = -S 3 y - X(S 3 yz + S 3 zy ) + 2 A 2 ( S\ 



Dim 



S: 



zzy) 



(5.41) 



As we can see, the two first terms in qi{X) are a direct generalization of the result for two 
sites, which hints at the existence of a general pattern. As for the N = 2 case we find once 
again, that even before attempting to solve (|3.18j) , we could have predicted from ( 3.20 ) 
that the matrices q r 2fc-i(A) can only be linear combinations of Sy,Sy Z ,S zy ,Sy yy ,S ZZ y and 
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S xxy (for k = 1, equation (|5.1|) tells us though that the coefficient of S xxy is zero. This will 
change for higher orders in perturbation theory). We can therefore write, 

q = a(X, k)S 3 + (3(X, k)(S 3 z + S 3 zy ) + 7(A, n)S 3 yyy + 6{\, n)S 3 xxy + e(A, K )S 3 zy , (5.42) 

where 

«( A > K ) = Yl7 i ^fc-iCA)^- 1 , /3(A, «) = 1 b 2k -i{X)K 2k -\ (5-43) 

^ — 'fc=l ^ — 'fc=l 

^ k ) = Er=i s 2fc-i( A )« 2fc_i . ^ «) = Er=i ^-i(a)^- 1 , (5.44) 

e(A, «) = J2" =i e 2fe _ 1 (A)« 2fe - 1 . (5.45) 
Computing coefficients up to order kJ we find the results in tables 5-7. 





-A u 


-A 


-A 


— A B 


-A 


-A iu 


-A 


ai(A) 


1 




















a 3 (A) 


i 

3 


8 
3 


16 














a 5 (A) 


1 
5 


122 
15 


144 


5024 
15 


3072 
5 








67(A) 


1 
7 


576 
35 


9616 
15 


432832 
105 


1755136 
105 


2720768 
105 


196608 
7 


rfi(A) 























d 3 (A) 








2 4 
3 














rfs(A) 





2 
3 


496 
15 


1184 
15 


2 w 
5 








d 7 {\) 





2 b 
15 


4432 
35 


86848 
105 


65024 
15 


754688 
105 


2 1B 

7 



Table 5: The coefficients a,2k+i (A) and d2fe+i(A) for k < 4. 





-A 


-A 


-A b 


-A' 


-A !J 


-A u 


-A 


MA) 

63(A) 
65(A) 
67(A) 


1 

4 
3 

23 
15 
176 
105 




28 

3 
664 

15 
4344 

35 




2 6 
3 

1568 
5 

13536 
7 





3328 
5 

52416 
5 





2 12 

5 

1104384 
35 






311296 
7 


0, 



0, 

2 18 

7 



Table 6: The coefficients 6 2 /c+i(A) for k < 4. 





-A y 


—A 


— A B 


-A" 


-A lu 


-A 


—A 


Sl (X) 


-2 




















53(A) 


-4 


-8 


2' 
3 














s 5 (A) 


_ 28 
5 


112 

5 


2592 
5 


4608 
5 


2 L3 

5 








S 7 (A) 


232 
35 


288 
35 


91008 
35 


452224 
35 


356352 
7 


491520 
7 


2 ia 
7 


ei(A) 


2 




















e 3 (A) 


20 
3 


56 
3 


2 V 

3 














e 5 (A) 


196 
15 


400 
3 


3872 
5 


6656 
5 


2 13 

5 








e 7 (A) 


440 
21 


53152 
105 


206336 
35 


75904 
3 


69632 


622592 
7 


2 iy 
7 



Table 7: The coefficients s 2 fc+i(A) and e2/c+i(A) for k < 4. 
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It is now possible to use these perturbative results to compute h(X, k) for particular 
values of A and n. We find that the structure of the Hermitian counterpart of the original 
Hamiltonian is: 

h(X, k) = i4 x (X, k)SI x + n 3 yy (X, K)Sy\ y + i4 z (X, n)Sl z + nl(X, n)Sl 

+/4r 2 ( A > K ) s Lz + /4/*( A > K ) S yyz + &zz(\ K )S 3 ZZZ , (5.46) 

which resembles the result for two sites, but includes few extra terms that couple all three 
sites. The functions n xx , . . . , [i zzz are all real functions of the couplings. As for N = 2, 
the Hamiltonian above is "PT-symmetric, which follows from the fact that all matrices 
involved are invariant under the adjoint action of the operator VT (see equation (|4.5| )). 
As for N = 2 also, these are the only matrices that are both VT symmetric and real 
(notice that, from the definition ( fP| ) for N = 3, it holds that S xxz = S zxx = S xzx and 

yyz zyy °yzy)' 

5.4 The N = 4 case 

It is interesting to investigate how the perturbative results generalize as we increase the 
number of sites. The N = 4 case is especially interesting as it is the simplest example 
for which we may see non local interaction terms in the Hermitian Hamiltonian. There is 



again only one solution for (71(A) which is compatible with the conditions ( 3.26 ), that is 



6A 3 lW 4 — <? 4 — <? 4 

u/x \>Jy UZ '-'yz zy 

'' ! ~ ' M "-" : " ' 40A 2 - 9 



<7i (A) — -S y - X(Sy Z + S zv ) 



+ 40A^9 [(9 " 32A2)A2(5 ^ + S ^y ] ~ 32X * S *y* ~ 2X2(3 ~ 16X ^ S tyy 

j4 1 c4 \ 1 o \ 3 1 c4 r c4 , c4 \ 



— 1\ 2 ( G I 4: — 2S" 4 4- S" 4 1 4- IX s ( S" 4 — ^S 14 4- < ? 4 

o/x \^> XX y 0,J xyx ' u yxx) ' " \ u xxyz ou xyxz ' u xxzy I 

+ 2X 3 {§Sy zzz - 7Sy yyz ) + 64X 5 (Sy yyz ~ S ■ ^ ' ^ 

In many ways, this is a simple generalization of the results of two and three sites. The 
matrices that enter the expression are to a large extent the same we find for less sites, but 
we have now extra contributions involving Pauli matrices sitting at all four sites of the 
chain, which was to be expected. There are however two major changes 

• the dependence on A of the coefficients is not polynomial anymore, 

• the first occurrence of non-local interactions appears through the matrix Sy UZ . 

As for lower values of N, it is not difficult to argue that the matrices ( |4.1[ ) entering 
the linear combination ( |5.47| ) are the only ones that are compatible with (3.24). Hence, 



as expected, the same structure extends to higher orders in perturbation theory, although 
expressions become extremely involved. The table below gives (73(A) as a sum of terms 
given by the matrices on the first column multiplied by the corresponding coefficients in 
the second column, 
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93 (A) 


Coefficients 


S 4 


-81+72A 2 + 1892 A 4 -4224A°+28672A 8 - 131072 A lu 
3(-9+40A 2 ) 2 


S 4 + S 4 

y z z y 


2916A-22842A 3 +27216A : '+81152A'+251904A'' + 786432A 11 -6291456A 13 
3(— 9+40A 2 ) 3 


S 4 

yuz 


64A 5 (351+276A 2 -8352A 4 -4096A 6 +65536A 8 ) 
3(-9+40A 2 ) 3 


C4 , o4 

y z z z zy 


4A 2 ( — 1215+7722 A 2 + 12432 A 4 — 15 1808 A 6 + 131072 A s — 262144A ln +2097152A 12 J 

3( — 9+40A 2 ) 3 


s 4 

zyz 


4A 2 (l215-16065A 2 +26952A 4 +72448A 6 -2097152A 12 ) 
3(— 9+40A 2 ) 3 


tit It 1 

yyy 


4A 2 (-729+8667A 2 -14040A 4 -97024A 6 +393216A s -1048576A lu +2097152A 12 ) 

3(— 9+40A 2 ) 3 


g4 , g4 

xxy yxx 


8A 2 (243-2079A 2 +7752A 4 + 16768A 8 -159744A S +262144A 1U ) 
3(-9+40A 2 ) 3 


o4 

xyx 


4A 2 (-972+6939A 2 -31512 A 4 +82176A 6 -188416A S +262144A 1U ) 
3( — 9+40A^) 


o4 , o4 

xxyz ' xxzy 


2A 3 (405-19368A 2 + 146048A 4 -349184A 6 - 131072 A B + 1048576A 1U ) 
3(-9+40A 2 ) 3 


o4 

^xyxz 


64A 3 (81-99A 2 -692A 4 +672A 6 +4096A 8 ) 
3(-9+40A 2 ) 3 


o4 

zyyy 


32A 3 (-567+6390A 2 -21448A 4 + 12096A 6 +49152A s -196608A 10 +524288A 12 ) 
3(-9+40A 2 ) 3 


a4 

yzzz 


32A 3 (-729+7290A 2 -19512A 4 -5952A e +32768A s -65536A ll) +524288A 12 ) 
3(-9+40A 2 ) 3 



In general we have, 



q = C(A, k)S* + 0{\, n){S% + S} z ) + #(A, k)S^ z 

+fx(X, k)(S^ zz + S A zzy ) + v(X, n)S% z + £(A, K)S^ yy + w(X, K)(S* xy + S^ xx ) (5.48) 



+g(X, K)S xyx + ?(A, K)(S xxyz + S xxzy ) + t(A, K)S xyxz + v(X, K)S zyyy + x(A, K)S yzzz , 
where all coefficients £(A, re), 6>(A, k), ... , x(A, k) can be expressed as expansions of the form 



( 5.43 ) and are real functions of the couplings. Perturbation theory results show that the 



Hermitian Hamiltonian h(X,K) has the following structure: 

h(X, k) = fj, xx (X, k)S xx + u xx (X, k)S xux + fiyy(X, K)S m + u yv (X, k)S, 



+A t zz(Aj K )S ZZ + v zz (X, k)S zuz + /x z (A, k)S z + fi xxz (X, k)(S xxz + S zxx ) 

^^■xzxi^i K )^xzx + Vyyzi^i K )(^yyz + ^zyy) + f^yzyi^i K )^yzy ^zzzi^t K )^zzz 
^l^xxxxi^i K )^xxxx ^yyyyi^i K ^yyyy ^zzzzi^i K )^zzzz ^xxyyi^i K )^xxyy 
^/^xyxyi^i K )^xyxy ~T" ^zzyy^i ^)^zzyy "T" /Azj/zj/ ( A > ^)^zyzy ~t~ ^xxzz^i K )^xxzz 
*xzxz \ 



+^xzxz(\ K ) S Lxz- ( 5 - 49 ) 



As expected from the expression of q, we find that h(X, k) involves non-local interaction 
terms proportional to S xux , S yuy and S zuz . The remaining terms are the natural general- 
ization of the those appearing for the iV = 2, 3 cases plus additional terms corresponding to 
interactions that couple all four sites of the chain. Once again, all coefficients [i xx , . . . , \x xzxz 
are real functions of the couplings. As for previous cases, it turns out that matrices ap- 



pearing in the linear combination ( 5.49| ) are exactly those that are both invariant under 



■pT-symmetry, according to equation ( |4.5| ), and real. 
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5.5 Some general features from perturbation theory 

We would like to end this section by summarizing the main results that we have obtained 
from our perturbative analysis. Since we have only solved for 2, 3 and 4 sites, our conclu- 
sions are based on a case-by-case analysis rather than rigorous proofs. However, we believe 
that the consistent occurrence of certain features across the various examples that we have 
studied provides strong support for these conclusions. 

Firstly we found that the combination of perturbation theory and the assumption of 
Hermiticity of the Dyson operator n = e q ^ 2 fix the metric p and therefore the Hermitian 
Hamiltonian h(X,n) with its corresponding observables uniquely. We have established this 
for N = 2, 3, 4 and arbitrary values of both coupling constants as well as for arbitrary N 
if A = 0. 

Secondly concerning the specific algebraic structure of the Hermitian Hamiltonian, we 
have seen that it becomes more involved for higher values of TV. For N > 2 it generally 
includes interaction terms that couple two or more adjacent sites, as well as non-local 
terms that couple non-adjacent sites. In addition, this structure is entirely dictated by 
VT symmetry, which selects out which tensor products of Pauli and identity matrices the 
Hamiltonian will be a linear combination of. Combining the requirement of VT symmetry 
with the requirement of h(X, k) being real completely fixes the general structure of h(X, k), 
although not the specific dependence on the coupling constants A and k, which is fixed by 
perturbation theory. All examples studied indicate that for a given value of N, all solutions 
(?2fc-i(A), with k > at different perturbative orders, share a common structure, namely 
they are all linear combinations of the same set of matrices, with coefficients that increase 
in complexity with increasing values of k. 

Finally, concerning the numerical accuracy of perturbation theory, we have demon- 
strated in detail that it converges very quickly for N = 2. For N = 2, 3 and 4 it becomes 
very difficult to perform computations up to such high orders of perturbation theory reached 
for N = 2 and the rate of convergence has not been analysed in detail for such cases. An 
interesting aspect of the model studied here is the dependence of the Hamiltonian on two 
coupling constants. For N = 2, we have carried out perturbation theory in both such 
couplings and found quick convergence in both cases. All our perturbation theory results, 
suggest that the entries of the Hermitian Hamiltonian h(X, k) can generally be expressed 
as a double Taylor series in A and k. 



6. Expectation values of local operators: form factors 

In this section we want to employ our general formulae in order to compute the expectation 
values of certain local operators of the chain. In particular, we will be looking at the 
expectation values of the total spin in the x and z directions in the ground state of the 
chain. These expectation values are commonly known as the magnetization in the x and z 
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directions. Recalling the results from section 3^, we define 

1 



M z (X,k) 
M x (X,k) 



(v g \vs?v\* 



(V g \vS?v\* 



Wg\S?\l/> g ) 
Wg\S*!\ll> g ) 



(6.1) 
(6.2) 



where \ip g ) is the ground state of the Hermitian Hamiltonian and |\P S ) is the ground state of 
the non-Hermitian one. We assume that the states are normalized to (i/jg\ip g ) = {^ g \^f g ) = 
1. In the following sections, we will carry out this computation for A = with generic N 
and for A 7^ for small values of N. 

6.1 General solutions for A = 



In section 4.1 we described in detail how for A = the original Hamiltonian and its 



Hermitian counterpart simplify greatly. Indeed, the latter can be found in all generality, 
for any number of sites, resulting in the expression ( [4 . 9| ) . Taking ( 4.10| ) and Q4.12 ) into 
account, it is very easy to show that 



N 

M z (0,k) = -, 



(6.3) 



which is nothing but the total spin of the chain and does not depend on the particular value 
of the coupling k. This result is to be expected for a Hamiltonian like ( [4.9| ). Naturally, the 
spins of the chain tend to align in the direction of the field, and will all be up so that the 
magnetization is just the total spin of the chain and maximal. A similar computation can 
be performed for M X (Q, k) for each particular value of N. In all cases one finds 



M x (0,k) = 0, 

which is also what one would expect for this model, as the Hamiltonian 
favour any particular direction of the spin a x . 

6.2 General solutions for k = 



(6.4) 



does not 



For k = the Hamiltonian ( |1.1[ ) is Hermitian and therefore computations of the magneti- 
zation simplify, as n = I. The ground state will nonetheless still depend on the value of A. 
For example, for N = 2 it is 



\i> 9 



2(1 + A 2 + x/l + A 2 ) 



/i + v / T+a i \ 



A 



(6.5) 



with energy E„ = — yl + A 2 and the magnetizations becomes simply 



M 2 (A,0) 



a/T+a 2 " 



and 



M a (A,0) = 0. 



(6.6) 
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The function M Z (X, 0) flows between the value M z (0, 0) = 1, as seen in the previous section, 
and M 2 (oo,0) — > 0. This is simply because for k = our model is nothing but the Ising 
chain with a magnetic field of intensity 1/A in the z-direction. Therefore, as A —* oo the 
intensity of the perturbing field tends to zero, and the ground state of the chain has zero 
magnetization, as consecutive spins align in opposite directions to minimize energy. This 
is a general feature that will also hold for higher values of N. For example, we find 

0) = i + — ; 2 ~ A = for N = 3, (6.7) 
2 2y /l + (-1 + A) A V ' 



1 _A 2 + v / l + A I Wl + A' + Vl + A 1 



M (A, 0) = ^ 7 for AT = 4. (6.8) 

\/2\A + A 4 

In both cases we recover the expected behaviour: /x(0, 0) = A 7 / 2 and //(oo,0) = 0. The 
K = curve in figure 4 is precisely a plot of the function fl6.6|) for N = 2. 

The second equation in ( |6.6j ) can also be explained easily as a consequence of the 
symmetry of the Hamiltonian H(X,0). Such Hamiltonian is invariant under the transfor- 
mation erf — > —af at each site i of the chain. This means that any form factor involving 
the operators erf must have the same symmetry. Therefore, 

M x (A,0) = -M,(A,0), (6.9) 

which implies M X (X, 0) = for all values of N. 

6.3 The N = 2 case for «, A ^ 

Let us now compute M Z (A, k) and M^A, k) in the more generic situation when both cou- 
pling A and k are non vanishing. We will start by analyzing the magnetization in the 
z-direction. In this case (A 7^ 0), the form of the ground state of the Hermitian chain 
is not particularly simple and therefore we will work with the first equality in ( |6.2[ ) and 
employ the properly normalized ground state of the non-Hermitian Hamiltonian. As figure 
3 shows, the magnetization is maximal at A = with value 1, and exhibits different kinds 
of behaviour as A increases, depending of the value of k under consideration. 

For every fixed value of k, the corresponding graph in figure 3 generally only covers 
a small region of values of A. These are precisely the values that lie in the region U-pj- of 
figure 1, namely those values for which all eigenvalues of H(X, k) are real. As shown in 
figure 3, the smaller the value of k the larger this region becomes in A. Depending on the 
value of k the magnetization exhibits a rich structure: for n > 0.7 it is a strictly decreasing 
function, whereas for k < 0.6 it has a minimum. This minimum is located near the critical 
value of A above which some eigenvalues of the Hamiltonian become complex, except for 
k = 0.6, where the minimum of the magnetization shifts to a smaller value of A. 

With regard to the magnetization in the x-direction we find that it vanishes for all 
values of A and k. This is so because the Hermitian counter-part of H(X, k) with N = 2 has 
the form ( 4.29j ) and therefore the Hamiltonian h(X,K) has the same symmetry described 



at the end of the previous section. 
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Figure 3: The magnetization in the z-direction for N = 2 as a function of A and k. 

It is also interesting to analyze how the presence of an imaginary magnetic field in 



the x-direction in (LI), as opposed to a real one really changes the physics of the model. 
Figure 4 precisely shows the magnitude of that change for the magnetization when k is an 
imaginary number. The Hamiltonian (1.1) is now that of the Ising spin chain with both 
a perpendicular and longitudinal fields applied at each site of the chain. The competition 
between these two fields will determine the values of the magnetization in both the x and 
z-directions. 




-10 -8 -6 




Figure 4: The magnetization in the x and z directions for N = 2 and k imaginary. 



We also observe that the magnetization is strictly smaller than 1, as it should be. 
Computing the expressions ( |6.1[) and (|6.2| ) in the standard metric p = I, i.e. disregarding 
the fact that the Hamiltonian is non-Hermitian, leads to non-physical values larger than 
one. 
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7. Conclusions 

We have demonstrated that there are various possibilities to implement 'PT-symmetry for 
quantum spin chains, either as a "macro-reflection" by reflection across the entire chain 
or as "micro-reflection" by reflecting at individual sites. These new possibilities constitute 
symmetries for the model H(X, k) in ( |1.1| ) we focussed on, i.e. Ising quantum spin chain 
in the presence of a magnetic field in the z-direction as well as a longitudinal imaginary 
field in the x-direction. However, there are also implications for other Hamiltonians such 
as Hxxz in Q2.10D and Hog in ( 2.12 ). Due to the various possibilities to implement parity 



the corresponding metric and therefore the underlying physical model is more ambiguous 
and it requires further clarification as to which physical system it describes. Remarkably 
the non-Hermitian Hamiltonian H(X, k) fixes the underlying physics uniquely under the 
sole assumption the Dyson map ry is Hermitian. As pointed out above this uniqueness is 
not obtained in general. One might conjecture that this is due to the finite dimensionality 



of the Hilbert space, as opposed to continuous models studied for instance in [16|, 29[, but 
our comments on Hxxz an d Hog suggest this is not the case. The explanation lies surely 
in the different types of symmetries a Hamiltonian might possess, which is supported by 
the fact that two different types of metric operators, say p and p, can always be used to 



define a new non- unitary symmetry operator S = pp 1 [34, 35 1. 

We have shown that all these possibilities serve to define anti-linear operators, which 
can not only be used to explain the reality of the spectra and identify the corresponding 
domains in the coupling constants, but can also be employed to define a consistent quan- 
tum mechanical framework. Regarding the technical feasibility of this programme, we have 
demonstrated for two sites that the perturbation theory, in k as well as in A, converges very 
fast by comparing it with the exact result. We took this as encouragement to tackle also 
three and four sites, albeit up to not as high orders of perturbation theory. Our pertur- 
bative analysis has allowed us to demonstrate for specific examples that the combination 
of perturbation theory and Hermiticity of the Dyson operator are sufficient to uniquely fix 
r], p and /i(A, k). In fact, for the model at hand, the constraint of Hermiticity of rj appears 
to be sufficient to entirely fix the algebraic structure of these quantities, even before any 
perturbative analysis is carried out. 

Clearly there are various open issues and follow up problems associated to our in- 
vestigations. Firstly one may try to complete the analysis for the Hamiltonian H(X, k) 
by carrying out further numerical studies, perturbative computations for more sites and 
ultimately obtain a complete analytic understanding for instance by means of the Bethe 
ansatz. Special attention should be given to the values of k and A corresponding to the 
exceptional points, when the usual analysis is expected to break down. Secondly one may 
consider the model for higher spin values as for instance studied in (7j. Finally it would 
be also very interesting to investigate some other members of the class belonging to the 
perturbed A^ Pi9 -series of minimal conformal field theories. 

Acknowledgments: A.F. is grateful to Giinther von Gehlen for bringing the papers H [7| 
to our attention. O.C.A. would like to thank Benjamin Doyon for helpful discussions and 
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A. Exact Hermitian Hamiltonian for N = 2 

As demonstrated in section |5| perturbation theory, both in k and A, agrees numerically very 
well with the exact results of section ||. We showed that the Hermitian counterpart to the 
Hamiltonian (|l . l| ) for N = 2 can be obtained by computing 



h(X,K) = e q/2 H(X,K)e' q/2 , with q = a(X, k)S 2 + (3(X, k){S 2 z + S; 



2 ) 

zyJ 



(A.l) 



where the functions a(A, k) and f3(X,n) have been evaluated perturbatively employing 
( -5.1 0|) in section || In terms of these functions and combinations thereof defined in ( 5.11| ) 
and ( [5 . 1 2| ) , the entries of the Hamiltonian (A.l) are given by, 



hn(X,k) 



h 2 2 (A, k) 

^44 (A, k) 



-(2 7 (A, k) 4 )" 1 \a(X, k)P(\, k) [6a(X, k)/3(A, «) - 7 (A, nf 
+2n sinh(7(A, n))a(X, k)j(X, n)5(X, «) 
+2 cosh( 7 (A, k))5(X, k) (2Aa(A, k)(3(X, k) + 6(X, «)) 
+Aa(A, k) (a(X, k) 2 — 3a(A, k)/3(\, k) + 4/3 (A, k) 2 ^ e(A, k 
+ cosh(27(A, k))q(A, k) ((3(X, re) — Aa(A, re)) e(A, re) 2 
+re sinh(27(A, re))/3(A, k)j(X, re)e(A, re) 2 
sinh(7(A, re)) 



7 (A,re) 



[sinh(7(A, re))a(A, re) (/3(A, re) — Aa(A, re)) 
+recosh(7(A, re))/3(A, re) 7 (A, re)] , 



(A.2) 



(A.3) 



(2 7 (A, re) 4 )" 1 a(A, re) ((6 - A)a(A, re)/?(A, re) 2 + (1 + 4A) /5(A, re) 3 

+a(A, re) 2 /3(A, re) (1 - 2A) - Aa(A, re) 3 ) + 2k sinh(7(A, re))a(A, re) 7 (A, re)<5(A, re) 
+2 cosh( 7 (A, re))<5(A, k) (2Aa(A, re)/3(A, re) + 6{X, re)) 
+ cosh(27(A, re))a(A, re) (Aq(A, re) — /3(A, re)) p(A, re) 2 



/ti 4 (A,re) 



— Ksinh(27(A, re))/3(A, K)"f(X, re)p(A, re) 
(2 7 (A, re) 4 )" 1 [-4cosh( 7 (A, re))a(A, re)/3(A, re) (2Aa(A, re)/3(A, re) + <5(A, re)) 

\2 



(A.4) 



+2re sinh(7(A, re))/3(A, re) 7 (A, re) ^— 2a(A, re) + cosh( 7 (A, re))<5(A, re) 
—p(X, re) ^Aa(A, re) 2 — 3a(A, re)/?(A, re) — 2A/3(A, re) 2 ^ e(A, re) 



cosh(27(A, re))a(A, re) (Aa(A, re) — (3(X, re)) e(A, n)p(X, re)] 



(A.5) 
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